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We discuss the effect of the thermal environment on the 
low-temperature response of the magnetization of uniaxial 
magnets to a time-dependent applied magnetic field. At 
very low temperatures the steps-wise magnetization curves 
observed in molecular magnets such as Mni2 and Fes display 
little temperature dependence where the apparent thermal 
assisted process are suppressed. However the changes of the 
magnetization at each step cannot be analyzed directly in the 
view point of a quantum mechanical nonadiabatic transition. 
In order to explain this deceptive apparent nonadiabatic be- 
havior, we study the quantum dynamics of the system weakly 
coupled to a thermal environment and propose a relation be- 
tween the observed magnetization steps and the quantum me- 
chanical transition probability due to the nonadiabatic tran- 
sition. 



PACS number: 75.40. Gb,76. 20. -Hq 



Magnetization processes of nanoscale molecules such 
as Mni2 and Feg have attracted much interest. For such 
small systems the discreteness of energy level plays an im- 
portant role and staircase structures of the response of 
the magnetization to a sweeping magnetic field have been 
observed The staircase is explained as a quantum 

mechanical transition at the avoided level crossing points, 
where levels of the Hamiltonian become almost degen- 
erate, and form repulsive structures as shown in Fig.l, 
which has been called resonant tunneling. This quantum 
mechanical transition has been studied from the point 
of view of the nonadiabatic transition |^ ll| . There are 
two characteristic features of each nonadiabatic transi- 
tion |lO| . One is the localization of the transition because 
it occurs only around avoided level crossing points. The 
other is the dependence of the transition probability on 
sweeping rate of the magnetic field, the energy gap and 
the gradients of the levels. Since at each avoided level 
crossing point only two levels play an important role, the 
transition probability can be described by well-known the 
Landau-Zener-Stiickelberg (LZS) mechanism [p^-[l4|]. 

However the behavior of these magnetic systems can 
easily be affected by thermal fluctuations even at low 
temperatures, because the energy scales involved are 
rather small. At relatively high temperatures (T ^ IK) 



the temperature dependence of the magnetization pro- 
cess is very significant, where excitations to higher levels 
provide other channels of resonance tunneling which is 
called thermally assisted resonant tunneling [T5|-[r^ . The 
external noise may affect the LZS mechanism itself which 
has been also studied JTsj-pOt. 

On the other hand, at very low temperatures (T ~ 60 
mK), the magnetization curve shows very little change 
with temperatures and only quantum mechanical phe- 
nomena seem to be dominant Q. However, as we will 
show below, even at such low temperatures thermal fluc- 
tuations cause inevitable effects which prevent a direct 
application of mechanism of the nonadiabatic transition. 

In this letter we investigate the effect of the thermal 
environment at very low temperatures on nonadiabatic 
transitions and find a relation between the observed data 
and the true quantum mechanical transition probability, 
from which the energy gap at the avoided level crossing 
point via the LZS formula can be deduced. 

Let us consider the change of magnetization when the 
external field is swept from a negative value to a posi- 
tive value. Initially the system is assumed to be in the 
ground state with the magnetization niQ ~ —S (approxi- 
mately). As the field increases, the state with mp crosses 
states with the magnetization S, S — 1, ■ ■ ■ , and 0. At 
each avoided level crossing point a nonadiabatic transi- 
tion occurs (Fig.l). 




Fig.l Schematic energy level diagram and the 
nonadiabatic transitions, pi denotes the probability 
that the system remains in the same eigenstate. 

We assign numbers i (i = 1, 2, ■ • ■) for the avoided level 
crossing point where the state of mo crosses a state with 
rrii ~ S—i + 1 (= 5*, 5* — 1, • • •, respectively). Let pi denote 
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the probability staying the same level at the ith avoided 
level crossing point. For pure quantum dynamical case, 
we have the following relation between the change of the 
observed magnetization at the crossing point i, AMi = 
Mi — Mi-i and the transition probabilities {pi}: 

i-l 

AM, = n (1 " ^'") " + " "^o} , (1) 

where Mi is the observed magnetization between avoided 
level crossing point i and i + 1. By this relation (|^), all 
the transition probabilities {pi} are obtained from the 
magnetizations in pure quantum cases. 

In the experiment of Pcrcnboom et al. for Mni2(5 = 
10) (T = 59mK) 0, shape of the magnetization process 
seems to saturate with the lowering of the temperature. 
When we analyze the data using the relation (|l|) , we can- 
not find any consistent set of the transition probabilities 
{pi}. In the experiment, the steps-wise changes of the 
magnetization occur at the avoided level crossing points 
where the state with mp ~ — 10 crosses with states with 
nii ~ 3, 2, 1 and {i = 7,8,9, and 10, respectively). 
The changes of the magnetization at the points are 0.62, 
3.54, 8.00 and 6.77, respectively. The relation (|l|) yields 
P7 = 0.0480, PS = 0.315, p9 = 1.13, and pw = -7.976, in 
contradiction to the trivial condition < Pi < 1. There- 
fore a naive application of nonadiabatic transition theory 
fails to explain the saturated magnetization curve in the 
very low temperature. 

We attribute this failure to the effect of thermal en- 
vironment even at such a low temperature. In terms of 
the potential picture (Fig. 2), the states with M > be- 
long to the right valley and we expect that these states 
easily relax to the bottom of the valley, i.e. to the state 
with M = S. Thus, once a quantum mechanical transi- 
tion from the metastable state of M = — S* to a state of 
M > takes place, the state is expected to relax easily to 
the lowest level due to some dissipation mechanism in the 
absence of an energy barrier. In the case of pure quan- 
tum transition, such a relaxation to the state of M = S 
is prohibited because the levels of the states are sepa- 
rated far away. If the time scale of dissipation is much 
shorter than that of the system and scale on which the 
magnetic field changes, the transfer to the lowest state 
takes a short time. As a result the magnetization curve 
will show a staircase as in the case of pure quantum dy- 
namics, but the change of magnetization at each step is 
difi'erent because of the relaxation transition M — > 5 in- 
stead of S" — i -I- 1. We call this steps- wise magnetization 
process in a dissipative environment "a deceptive appar- 
ent nonadiabatic magnetization process" . 



M=-S 



M=S 



Fig. 2 Potential picture of the metastability. 

In this scenario, we assume the following three prop- 
erties: (i) First quantum mechanical transition for mo(~ 
— S*) — > mi occurs with the probability of the pure nona- 
diabatic (LZS) transition {p]"^^}, and then (ii) the relax- 
ation from TTii — !■ mi(~ S) occurs by some dissipation 
mechanism, (iii) There is no relaxation directly from toq 
by the dissipation mechanism and thus the amount of 
magnetization change depends only on {pf^^} and does 
not depend on the temperature. Replacing rrii by toi in 
the relation (|l|), the change of the magnetization in this 
case is given by 

AMi = 11(1 -p„){[mo(l - pi) +mipi] - Too}- (2) 

n=l 

Using the data of now yields a reasonable solution for 
the {p^ys■. P7 = 0.0313, P8 = 0.185, pg = 0.515, and 
Pw = 0.898. 

In order to demonstrate that the above three prop- 
erties are really possible at very low temperatures, we 
simulate a relaxation phenomena of a magnetic system 
which very weakly couples to the external bath. Here we 
use a quantum master equation [ pT[ , 



dp{t) 
dt 

where 



I [H, p{t)] - A ([X, Rp{t)] + [X, Rp{t)]^) , (3) 



(fc|i?|TO) = Ci^^-J^MEk - i?„)(fc|X|TO), 

(^{lu) — I{uj) — I{~uj), and n^(tj) 



Here p is an inverse temperature of the reservoir l/T, 
and we set ?i = 1. \k) and |to) are the eigenstates of H 
with the eigenenergies Ek and Em, respectively. I{uj) is 
the spectral density of the boson bath. We take here the 
infinite number of phonons with the Ohmic dissipation 
I{uj) Iqlu |2^. As a more realistic bath for the exper- 
imental situation at very low temperature, we may take 
the dipole-field from the nuclear spins or other types 
of spectrum such as super-Ohmic type. X is an opera- 
tor of the magnetic system that interacts linearly with 
the bosons of the reservoir. The relaxation process can 
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be affected by the form of interaction of the system with 
the thermal bath, i.e. by the choice of X . Here we take 
X = \{S^ + S^). Generally X = S^^ is more efficient 
than X — Sz for the relaxation. A detailed comparison 
with other choices will be presented elsewhere. The alter- 
nate choices of concrete form of thermal bath, however, 
do not cause any significant qualitative change because 
the couplings to the bath is very weak. 

For Mni2, detailed form of the Hamiltonian has been 
proposed However the energy gap of the Mni2 is 

too small to observe the phenomena within the available 
computation time. Thus, here, we demonstrate the qual- 
itative features of the dynamics, i.e. the three properties 
(i), (ii), and (iii). We believe that the key ingredients 
of the general qualitative feature are the existence of the 
avoided level crossing points and weak coupling to the 
external bath. For the realistic model with much small 
energy gap, the features observed here should be realized 
in a much longer time scale. Thus we adopt a minimal 
model of a uniaxial = 10 spin system with the two 
ingredients: 



(4) 



with a linearly increasing external field, i?ext = ct — Hq 
where c is the sweeping velocity. The transverse field F 
represents the terms causing quantum fluctuations. We 
choose D = 0.1, F = 0.5 throughout this letter. 

In order to see the difference of relaxations between the 
case with and without the potential barrier, we compare 
typical two cases: (1) ffcxt — 0.05 and (2) i/cxt = 0.15 
and set the sweep velocity c — 0. As the initial state we 
take the second level, as indicated in Fig. 3. 
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Fig. 3 Energy level diagram of the model (Q) as a 
function of i/cxt- The white and black diamonds 
correspond to the case (1) and the case (2), respectively. 

The second level has M ~ —10 in the case (1) and M ~ 9 
in the case 2. In the both cases, the ground state has 
M ~ 10. The parameters are set to T = 0.1, Jo = 1.0, 
and A = 1.0 X 10"''. We study the relaxation for both 
cases by solving Eq. (||) . These probabilities are given by 
a diagonal element of p{t), i.e., (l|p(t)|l) and (2|p(t)|2), 
respectively. We observe almost no damping in the case 
(1), whereas a rather fast relaxation occurs in the case 



(2). Thus at a fairly low temperature, the thermal en- 
vironment causes significant difference in the relaxation 
process depending on the presence of a potential barrier. 
The difference between the cases (1) and (2) can be un- 
derstood analyzing the matrix elements of Eq. (^) . 

We now investigate time evolution of the system for a 
sweeping field c = 1.0 x 10~^ starting at Hq — —0.05. 
We study the case of pure quantum dynamics (A — 0)[P] 
and the case with a weak dissipation (A = 1.0 x lO""*) 
[D]. These magnetization curves are shown in Fig. 4. We 
show data for iJcxt > 0.45 because almost no change is 
observed for iJcxt < 0.5. 
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Fig. 4 Magnetization as a function of iJoxt- The dashed 
line denotes the pure quantum dynamics [P], and the 
solid line denotes the dissipative quantum dynamics [D]. 

For the case [P], we observe oscillating behavior due to 
spin precession, whereas in the case [D] this detailed 
structure is smoothed out by the dissipation. We find 
steps-wise magnetization curves in both cases. The 
changes of the magnetization are listed in Table I. 



Cross point (mo, rrii) 






(-10,5) 


0.511 


0.693 


(-10,4) 


8.32 


13.3 


(-10,3) 


3.50 


5.03 



Table I The changes of magnetization. AAf[p] and 
AM[p] are the changes for the case [P] and [D], 
respectively. 
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{mo, mi) 


P\R],i 


P\p].i 


P\'D],i 


P\D].i 


P\'D],i 


6 


(-10,5) 


0.0280 


0.0341 


0.0460 


0.0346 


0.0291 


7 


(-10,4) 


0.730 


0.616 


0.995 


0.688 


0.716 


8 


(-10,3) 


1.000 


0.726 


78.0 


0.835 


0.970 



Table II The transition probabilities obtained by 
various ways, see the text. 

From these data we estimate the transition probabili- 
ties by the relation (|l]) and (|^) which are listed in Table 
II. First we obtain the transition probabilities from the 
data in Table I setting too = —S and to^ = S—i+1. From 
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the data AM[£)], unacceptable probabilities {p[D],i} are 
deduced by the relation (|l|), while acceptable ones {p[D],i} 
are obtained by the relation (0). {p[D],i} agree with 
{P[p],i} obtained by the relation from the data AAf[p]. 
This agreement shows the three properties (i), (ii), and 
(iii) really realized in the present model and thus we can 
estimate the quantum mechanical nonadiabatic transi- 
tion by the relation (||). Although the magnetization 
TTii is almost constant: mg ~ — S*, ~ S* — i + 1 (i > 1), 
they show a little dependence on the magnetic field i?ext- 
Taking the iJext dependence of rrii into consideration, we 
also calculated the transition probabilities in the case [D] 
with (0). They are shown as {p[D],i}- We confirmed that 
{P[u],i} agree with the probabilities {pp,],^} directly ob- 
tained from the diagonal elements of the density matrix. 
The difference between ppj j andppj j simply come from 
the large value of T for the convenience of simulation. If 
r is very small as the case of the experiment, mi is very 
close to S — i + 1 and it is expected that P[D],i and P[D],i 
are very close. We present the time evolution of {i\p\i) in 
Fig. 5. This figure explicitly shows the three properties 
(i), (ii), and (iii). 

1.2 I • • • . • • • ^ 
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Fig. 5 The time evolution of the probability of individual 

levels. 

We estimate the energy gap from the transition prob- 
abilities with the extended LZS formula p^^: 



1 — exp 



7r(AS, 



2{mi — mo)c 



(5) 



where c is the changing rate of the Zeeman energy. Using 
{P[D],i}i we obtain the energy gaps for the avoided level 
crossings as AEq = 1.83 x IQ-^, AEr = 10.1 x IQ-^. 
These estimates agree with the correct value AEq = 
1.54 X 10"3 and AE^ = 10.0 x 10"^ directly obtained 
from the energy levels [|5|. If we use ppj.ij we have, 
of course, almost complete values, AEq = 1.57 x 10"^, 
AEr = 9.9 X 10~^. Thus we conclude that we can esti- 
mate the energy gap from the deceptive apparent mag- 
netization by the relation (^. 

In summary, we have considered a mechanism for 
deceptive apparent nonadiabatic magnetization process 



which is relevant when the temperature is very low and 
no temperature dependence is observed apparently and 
propose the general relation (^) between the steps in the 
magnetization and the energy-level splittings at very low 
temperature. With the relation we have estimated the 
quantum transition rate {pi} in the low temperatures 
. We demonstrated an example of deceptive apparent 
nonadiabatic magnetization process in a minimal model 
with the avoided level crossing points and weak coupling 
to the external bath. Elsewhere we will report on our 
investigation of the energy gaps {AEi} of Mni2 and Feg 
based on the detailed information of the values of jumps 
and the scanning speed c. 

The present study is partially supported by the Grant- 
in Aid for Research from the Ministry of Education, Sci- 
ence and Culture. 



[1] 
[2] 
[3] 

[4; 

[5 

[6: 

[7 

[8 
[9 

[lo; 

[11 

[12 
[13 
[14 
[15 

[16 

[17 

[18 

[19 

[20 
[21 



J. R. Friedman and M. P. Saracliik, T. Tejada and R. 

Ziolo, Pliys. Rev. Lett. 76, 3830 (1996). 

L. Thomas, F. Lionti, R. Ballou, D. Gatteschi, R. Sessoli 

and B. Barbara, Nature 383, 145 (1996). 

J. M. Hernandez, X. X. Zhang, F. Luis, and T. Tejada, 

J. R. Friedman, M. P. Sarachik and R. Ziolo, Phys. Rev. 

B 55, 5858 (1997). 

L. Thomas et al.. Nature 383, 145 (1996). 
F. Lionti, L. Thomas, R. Ballou, Barbara, A. Sulpice, R. 
Sessoli and Gatteschi, J. Appl. Phys. 81, 4608 (1997). 
C. Sangregorio, T. Ohm, C. Paulsen, R. Sessoli, D. Gat- 
teschi, Phys. Rev. Lett. 78, 4645 (1997). 
J. A. A. J. Perenboom, J. S. Brooks, and S. Hill, and T. 
Hathaway, and N. S. Dalai, Phys. Rev. B 330 (1998). 
S. Miyashita, J. Phys. Soc. Jpn. 64, 3207 (1995). 
S. Miyashita, J. Phys. Soc. Jpn. 65, 2734 (1996). 
H. De Raedt, S. Miyashita, K. Saito, D. Garcia-Pablos, 
and N. Garcia, Phys. Rev. B, 56 11761 (1997). 
S. Miyashita, K. Saito, and H. De Raedt, Phys. Rev. 
Lett. 80, 1525 (1998). 

L. Landau, Phys. Z. Sowjetunion 2, 46 (1932). 

C. Zener, Proc. R. Soc. London, Ser. A137, 696 (1932). 

E. C. G. Stiickelberg, Helv. Phys. Acta 5, 369 (1932). 

D. A. Granin and E. M. Chudnovsky, Phys. Rev. B 56, 
11102 (1997) 

A. Fort, A. Rettori, J. Villain, D. Gatteschi, and R. Ses- 
soU, Phys. Rev. Lett. 80, 612 (1998) 

F. Luis, J. Bartolome, and F. Fernandez, Phys. Rev. B 
57, 505 (1998) 

Y. Kayanuma and H. Nakayama, Phys. Rev. B 57, 13099 
(1998) 

V. V. Dobrovitski and A. K. Zvezdin, Europhys. Lett. 
38, 377 (1997) 

L. Gunther, Europhys. Lett. 39, 1 (1997) 
K. Saito, S. Takesue, S. Miyashita, unpublished, cond 



mat/9810069 



[22] H, Grabert, P. Schramn, and G. Ingold, Phys. Rep. 3, 



4 



115 (1988). 

[23] N. V. Prokof 'ev and P. C. E. Stamp, J. Low Temp. Phys. 

104, 143 (1996). 
[24] M. I. Katsnelson, V . V. Dobrovitski and B. N. Harmon, 

|cond-mat/9807174 
[25] Due to the large value of T and thus pg ~ 1, it is difficult 

to estimate precisely the concrete value of AEs- 



5 



